Unitarity-Limited Elastic Collision Rate in a Harmonically- Trapped Fermi Gas 
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We derive the elastic collision rate for a harmonically-trapped Fermi gas in the extreme unitarity 
limit where the s-wave scattering cross section is o'(fe) = 47r/fc^, with hk the relative momentum. 
The collision rate is given in the form F = 'yI{T/TF) — the product of a universal collision rate 
7 = ksTp / {Qirh) and a dimensionless function of the ratio of the temperature T to the Fermi 
temperature Tp. We find I has a peak value of ~ 4.6 at T/Tp ~ 0.4, 7 ~ 82 {T/Tpf for T/Tp < 
0.15, and I ~ 2{Tf/T)^ for T/Tp > 1.5. We estimate the collision rate for recent experiments on a 
strongly-interacting degenerate Fermi gas of atoms. 

PACS numbers: 03.75.Ss, 32.80.Pj 
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INTRODUCTION 

Recently, we and several other groups |M IE ^ Hi 
have begun exploring the strongly-interacting regime in 
degenerate Fermi gases of atoms. In these experiments, 
a magnetic field is applied to the atomic samples to 
tune the interparticle interactions to the vicinity of a 
Feshbach resonance where the scattering length is large 
compared to the interparticle spacing. In this regime, 
new forms of high temperature superfluidity are pre- 
dicted d, Qj and strongly anisotropic expansion has 
been observed 0, 0, . As we pointed out in Ref. , 
the strongly interacting regime leads to unitarity-limited 
mean field interactions as well as unitarity-limited colli- 
sion dynamics. In the latter case, the scattering cross sec- 
tion is the order of in/kp, where fci;- is the Fermi wavevec- 
tor. In the unitarity limit, the collision rate assumes a 
universal form and is proportional to the Fermi energy 
ksTF- At sufficiently low temperatures, Pauli blocking 
may suppress the unitarity-limited elastic collision rate 
for the trapped gas, producing an effectively coUisionless 
regime. However, a theoretical study of Pauli blocking 
in the unitarity-limited regime has not been presented 
previously, making it difficult to accurately estimate the 
collision rate. The primary purpose of this paper is to 
present such a treatment. 

Several groups have examined the effects of Pauli 
blocking on the elastic collision rate for an energy- 
independent cross section 0, . For comparison to 
the collision rates obtained with a unitarity-limited cross 
section, we begin by deriving a formula for the collision 
rate in a harmonic trap as a function of temperature 
for an energy-independent cross section a. The results 
reproduce those obtained in Ref. within 10% p^ . 
We then extend the treatment to include the energy de- 
pendence of the cross section in the extreme unitarity 
limit, where the zero energy scattering length 05 satis- 
fies |A:i;-as| ^ 1 and the gas is strongly interacting. We 
show that the numerically calculated collision rates for 
both the energy-independent and unitarity-limited cross- 
sections agree with analytic expressions derived for the 



high-temperature limit. 

The final part of the paper defines a hydrodynamic pa- 
rameter (f) = T/luj_, the ratio of unitarity-limited collision 
rate F to the transverse oscillation frequency u!±_ of atoms 
in the trap. We calibrate <j) by observing the threshold for 
hydrodynamic expansion of a strongly-interacting Fermi 
gas as a function of evaporation time |l|. We estimate 
(j) for several recent experiments on strongly interacting 
Fermi gases 0, 0, |B| . 



CALCULATING THE COLLISION RATE 

We consider a sirnple model that assumes s-wave scat- 
tering is dominant [l3|. In this case, the collision cross 
section takes the form 



a(fc) = 



Ana's 



(1) 



where k is the relative wavevector of a colliding pair of 
spin up and spin down fermionic atoms. 

In the trap, the average collision rate per particle, F, 
is determined from the s-wavc Boltzmann equation |l4| 
under the assumption of sufficient ergodicity. We con- 
sider the rate for the process in which a spin up and a 
spin down atom of total energy ein — £3 + ei collide to 
produce atoms with total energy eout = ei + £2- The ef- 
fects of Pauli blocking are included for the particles on 
the outgoing channel, and we assume a 50-50 mixture of 
atoms in the two spin states. The depletion term in the 
Boltzmann equation for the particle of energy 64 is inte- 
grated over £4 to determine the collision rate F for either 
spin state (as a collision inherently includes one atom of 
each spin). For an energy- independent cross-section, the 
integrated loss rate is then N/2 = —TN/2, and 

N Ma f 
^TT = ^773- / deide2de3de4'D{€min) x (2) 

S{ei + £2 - £3 - e4)(l - /l)(l - /2)/3/4, 

where F is the number of collisions per second per atom 
and N is the total number of atoms in the trap. Here, 
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I^(einin) is the density of states evaluated at the energy 
emin = min{ei,e2,e3,e4}- fi = ^/{gi + 1) is the occu- 
pation number with q,; = exp[(ei — ii)/kBT], and /i the 
chemical potential 15]. At zero temperature, the chem- 
ical potential is given by the Fermi energy fi(0) = ep = 
{SNy^^hu! = ksTp, where ui = {uj^uj^Y/^ with and 
Wz the transverse and axial trap oscillation frequencies of 
a cylindrically symmetric trap. 

For fermions, the collision cross section of Eq. ^ for 
\kas\ <C 1 is half that for bosons, i.e., a — 47ra|. We 
begin by determining F for this case. 

Energy-Independent Cross-Section 

The integrand in Eq. |2] is readily shown to be sym- 
metric under the interchange of all four particle labels. 
Hence, without loss of generality, we multiply the in- 
tegrand by 4, and take emin — ci, and I?(emin) — 
[e\/{21%^u)'^y\e2i9j,ieii, where 6*21 = 9{e2 - ei) is a unit 
step function. 

It is useful to write the collision rate as the product of a 
natural collision rate 7E1, which depends on the trap pa- 
rameters, and a dimensionless integral I-Ei(T/Tp), which 
describes the temperature-dependence. 
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FIG. 1: Temperature dependence of the elastic scattering col- 
lision rate F in units of the natural collision rate 7. The 
dashed lines indicate the high- and low-temperature approx- 
imations. A) The collision rate for an energy-independent 
cross-section. B) The collision rate for a unitarity-limited 
cross-section. 



r = 7Ei^Ei(r/rF) 



(3) 



We take the natural collision rate to be the classical 
collision rate at T — Tp, 



7EI 



NMoj^a 

ilT^kBTp' 



(4) 



Note that the rate is 1/4 of that obtained in a spin- 
polarized Bose gas. With this choice, /ei becomes 



Unitarity-Limited Cross-Section 

To include the energy dependence of the cross-section, 
we adopt the notation of Rcf. 14], and make the replace- 
ment 

27rM f /-f+W 
gP(£min) / dx dP(j{q), (6) 

[ZnTi)'' Jt/(x)<£„i„ Jp-{^) 



Im{T/Tp) - 144 



dxi dx2 dx3 x\ f{xi + X2) ^ 



f{xi + Xs) [1 - /(Xl)] [1 - /(Xl +X2+ X3)]. (5) 

Here f{x) = l/[g{x) + 1], where g{x) = exp[(TF/T)(x - 
IJ^/ep)]. We assume that for the cases of interest, the trap 
depth is large compared to ep and ksT. 

/ei is readily determined by numerical integration us- 
ing standard results for the chemical potential as a func- 
tion of T/Tp At low temperature, T/Tp < 0.2, 
we find that /ei is well fit by I-Ei{T/Tp) ~ 15 {T/Tpf, 
which displays the quadratic dependence expected for 
Pauli blocking in both final states. At high tempera- 
ture, T/Tp > 1.5, we find the expected temperature- 
dependence, Iei{T/Tp) ~ Tp/T as shown below. The 
complete function ImiT/Tp) is plotted in Fig. UK. The 
maximum value, /ei — 1.3 occurs for T/Tp ~ 0.5. This 
demonstrates that 7E1 is essentially the maximum colli- 
sion rate. 



where 2q = xl Pi- + Pi — P"^ determines the rela- 



tive wavevector q, P± = y/2M[e — e,nin — U{x.)] ± 
Y/2A/[ei„in — U{x)], and e is the total energy of the col- 
liding particles. 

We are interested in the extreme unitarity limit, where 
a{k) = Aw/k"^ according to Eq. Qand the elastic collision 
rate is the maximum possible. We write the collision rate 
as 



r = 7UL IvhiT/Tp). 



(7) 



The natural collision rate, Eq. 2]with cr = Air/kp, where 
= 2MkBTp/h^, then takes the fo rm 



7UL 



(8) 



The dimensionless integral /ul is similar to that of 
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Eq.m 

nOO nOO nOC 

Ivl{T/Tf)^IU / dxi / dx2 / dxsxjx 
Jo Jo Jo 

f{xi + X2) fixi + X3) [1 - /(a;i)] X 

[1 - f{xi +X2+ X3)] F{2xi +X2+ X3, xi), (9) 

where F{x,Xm) determines the energy-dependent cross 
section a{q) = Air/q^ in units ofAn/kp. The arguments 
of F{x,Xm) are x — e/tp-, and x,„ = e-aiin/^F- As in 
Eq. [SJ we take Xm = xi without loss of generahty. 
For a harmonic potential, 
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du u 



■n^lxra Jo Vx~~2 



In 



X- ia,u) 



where x±{'^i ^) — ^(cHj ''A ^ ^(Q!, m), ^(«, u) = a + 2(1 — 
2u^), B{a,u) — 2^/2(1 — u'^){a — 2u'^), and a = x/xm- 

We can simplify the form of the integrals in Eq. |51 by 
transforming from coordinates {xi,X2,X3,u} into coor- 
dinates {w,y,z,s}, where w = xi, y = {x2 + X3)/2xi, 
z = (x3 — X2)/2xi, and s — \ ~ v?'. After the transfor- 
mation we have 

poo poo 

Ivl{T/Tf)^1U dyG{y) dww^l-f{w)]x 
Jo Jo 

[1 - f[w{l + 2y)]] r dz f[w{l + y + z)]x 
J-y 

f[w{l + y~z)] (10) 

with G{y) given by 



G{y) = - f ds 

Jo 



1 - s 



Ms + y/2 + ^s{s + y)) - \n{y/2)]. (11) 

Eq. ^|is plotted in Fig. ^3. The maximum collision 
rate is larger by a factor ~ 4 than for an energy inde- 
pendent cross section a — An/kp, consistent with esti- 
mates from radio-frequency measurements of the mean 
field shift in Ref. For T/Tp < 0.15, we again 

find the quadratic temperature-dependence that results 
from Pauh blocking, /ul ^ 82{T/Tf)'^. In the high- 
T/Tp limit we find Jul ^ 2 {Tp/Tf. This matches the 
high-temperature prediction given below. The maximum 
value of Jul — 4.6 occurs at T/Tp ~ 0.4. 



and Pauli-blocking can be neglected. Including the de- 
pendence of the scattering cross section on the relative 
speed Vj., we find generally that 

N f 

rHTY= / o?xn|(x)n|(x)(u,.cr(u,.)), (12) 

where the angled brackets denote an average over the 
relative velocity distribution for pairs of atoms and 
J o?xn|,|(x) = N/2. Eq. E|has a simple physical inter- 
pretation: The spin up atoms at position x, i.e., n|(x) c?x, 
are hit by spin down atoms at a rate ni(yi){vrcr(vr)) ■ For 
an energy-independent cross section, we obtain Fht = 
7ei Tp /T in agreement with the numerical results for 
Iei{T/Tp) of Fig. The temperature dependence in 
this limit arises from the flux, which is the product of the 
density and the relative velocity, i.e., n(u,.) cx 1/T. For 
the extreme unitarity-limited cross section a = 47r/fc^, 
with = hk, we obtain Fht = 7ul 2(7^/2^)^ in agree- 
ment with the numerical results for IvhiT/Tp) shown in 
Fig. ^3. In this case, the 1/T^ temperature dependence 
arises because both the flux and the cross section vary as 
1/T. 



HYDRODYNAMIC PARAMETER (jj 

The collisional state of a gas can be described by a 
hydrodynamic parameter <j> which is the number of col- 
lisions that an atom experiences during a characteristic 
timescale. When cj) is large, the gas is coUisionally hydro- 
dynamic, while when cj) is small, the gas is collisionless. 
The actual values of cj) which qualify as "large" or "small" 
are determined by calibration. 

We choose the characteristic time-scale to be l/oj, 
and take (jj — T/uj where lu — lu_i_,u!z are the oscilla- 
tion frequencies of atoms in the trap. These are also 
the natural time scales for ballistic expansion, where the 
size of the cloud scales as ^/T+l^J^ in each direction. 
In the unitarity-limited regime, and for a cylindrically- 
symmetric trap with elongation parameter A = w^/wj,, 
we can write (f)± as 



')± = 7 IvL[T/Tp) 

DTT 



(13) 



with iV the total number of atoms in the 50-50 mixture. 
Then, (/)^ = 



Comparison With Analytic High-r/Tp Results 

We check the numerical results for the temperature 
dependence of the rates by calculating the collision rate 
Fht for T ^ Tp directly from the phase-space s-wave 
Boltzmann equation [li | in the high-T/Tp hmit, where 
the occupation number is given by a Boltzmann factor 



Calibrating 

We now turn to the question of determining the ap- 
proximate value of 4> for which the transition between 
collisionless and collisional behavior occurs. In Ref. 0], 
we investigated the anisotropic expansion properties of a 
strongly-interacting, degenerate Fermi gas of ^Li. When 
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FIG. 2: Observed aspect ratio of a strongly-interacting Fermi 
gas after 600 /is free expansion as a function of evaporation 
time. The aspect ratios corresponding to ballistic and hydro- 
dynamic expansion are indicated by the dashed, horizontal 
Unes. The solid curve has been included to guide the eye. 
The gas evolves smoothly from ballistic expansion to hydro- 
dynamic expansion as the evaporation time is increased. 

released from a highly-elongated trap, the originally nar- 
row dimensions of the gas expanded rapidly, while the 
broad dimension remained largely unchanged — inverting 
the aspect ratio of the cloud. We have studied how the 
observed aspect ratio of the expanded cloud varies with 
the duration of evaporative cooling. The aspect ratios 
are measured for a fixed expansion time of 600 /iS and 
are compared to the predictions of ballistic and hydro- 
dynamic expansion. The results are plotted in Fig. |21 

From this figure, we see that the hottest clouds (short 
evaporation times) expand ballistically, while the coldest 
clouds (longest evaporation times) expand hydrodynam- 
ically. Ballistic expansion is expected in a normal, coUi- 
sionless gas. Since the rapid transverse expansion extin- 
guishes collisions before the axial distribution can change 
significantly, the coUisional behavior of the expanding gas 
can be associated with (px . For the shortest evaporation 
times shown in Fig. |21 we observe ballistic scaling with 
T = 3Tf and = 4 X 10^. For our trap A = 0.035. From 
Eq. 1131 we obtain 0_l — 0.4 which then corresponds to 
approximately collisionless behavior. Therefore, (j) — 0.4 
is, in general, the condition for collisionless behavior for 
any timescale l/iv. 

Application to Experiment 

In addition to our work in , several groups have also 
recently observed hydrodynamic expansion of a Fermi gas 
in the strongly-interacting regime 0, . The authors of 
these papers claim that their experiments are collisionally 
hydrodynamic. For Ref. we take the total number of 
atoms iV = 2.4 X 10^ and T/Tp = 0.34, with A = 0.016, 
yielding Iul = 4.5 and 0x = 5.4. For Ref. [3, we take 
iV = 7 X lO'^, T/Tf = 0.6 and A = 0.35, yielding Iul = 



3.5 and (f>±^ = 7.7. Hence, we agree with their conclusions. 

Our experiments (l| produce strongly-interacting 
Fermi gases at considerably lower temperatures. In those 
experiments, the total number of atoms is ~ 1.5 x 10^ 
for temperatures 0.08 < T/Tp < 0.2. Eq. US yields 
0.7 < 0_L < 3.7, while 20 < 0^ < 106 . The values of 
corresponding to our lowest temperatures indicate that 
the trapped gas is nearly collisionless on the transverse 
timescale, but collisional on the axial. 

The onset of high-temperature superfluidity has been 
recently predicted in the temperature range T/Tp = 
0.25 — 0.5 [1,0, 111- Since Pauli blocking is ineffective for 
the unitarity limited cross section when T/Tp > 0.25, it 
is not clear how collisions in the normal component will 
affect the formation of this high temperature superfluid. 

For an expanding gas, we cannot make a definitive 
statement about the collisional nature, even if it were col- 
lisionless when trapped, as Pauli blocking may become 
ineffective as a result of nonadiabaticity in the expan- 
sion, deformation of the Fermi surface, or through other 
effects 01 . The magnitude of these effects, and to what 
extent they modify (j) remains an open question. We are 
therefore working on experiments which will directly de- 
termine if the gas contains a superfluid fraction. 
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